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Abstract. 

Any two-dimensional infinite regular lattice G can be produced by tiling the plane 
with a finite subgraph B C G; we call B a basis of G. We introduce a two-parameter 
graph polynomial Pb (<?, v) that depends on B and its embedding in G. The algebraic 
curve Pe(g,w) = is shown to provide an approximation to the critical manifold of 
the g-state Potts model, with coupling v = e K — 1, defined on G. This curve predicts 
the phase diagram not only in the physical ferromagnetic regime (v > 0), but also 
in the antiferromagnetic (v < 0) region, where analytical results are often difficult to 
obtain. For larger bases B the approximations become increasingly accurate, and we 
conjecture that Pb(q, v) — provides the exact critical manifold in the limit of infinite 
B. Furthermore, for some lattices G — or for the Ising model (q — 2) on any G — the 
polynomial Ps(g, v) factorises for any choice of B: the zero set of the recurrent factor 
then provides the exact critical manifold. In this sense, the computation of Ps{q,v) 
can be used to detect exact solvability of the Potts model on G. 

We illustrate the method for two choices of G: the square lattice, where the 
Potts model has been exactly solved, and the kagome lattice, where it has not. For 
the square lattice we correctly reproduce the known phase diagram, including the 
antiferromagnetic transition and the singularities in the Berker-Kadanoff phase at 
certain Beraha numbers. For the kagome lattice, taking the smallest basis with six 
edges we recover a well-known (but now refuted) conjecture of F.Y. Wu. Larger bases 
provide successive improvements on this formula, giving a natural extension of Wu's 
approach. We perform large-scale numerical computations for comparison and find 
excellent agreement with the polynomial predictions. For v > the accuracy of the 
predicted critical coupling v c is of the order 10~ 4 or 10~ 5 for the 6-edge basis, and 
improves to 1CT 6 or 10 -7 for the largest basis studied (with 36 edges). 



We dedicate this article to Professor Fa-Yueh Wu on the occasion of his 80th birthday. 
1. Introduction 

The two-dimensional g-state Potts model [1] has been extensively investigated over 
the last sixty years and has served as a testbed for the powerful tools of quantum 
integrability [3J. It is defined by the reduced hamiltonian 
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where K is the dimensionless coupling between the spins cxj = 1,2, ... ,q that are defined 
on the the vertices i G V of some given connected graph G = (V, E) with vertex set V 
and edge set E. Here (3 = l/k^T, where T is the temperature and k-Q the Boltzmann 
constant, and £(crj, cr,) is Kronecker's delta function. From now on we suppose that G 
is a piece of a regular lattice, and we are primarily interested in the limit of an infinite 
two-dimensional lattice. 

The partition function Z corresponding to (Tj[|) is most conveniently expressed in 
the Fortuin-Kasteleyn representation [lj 

z = ^M-m) = $>iv (A) , (2) 

a ACE 

where \A\ denotes the number of edges in the subset A, and k(A) is the number of 
connected components in the induced graph Ga = (V,A). The temperature parameter 
v — e K — 1 will be used instead of K in the following. Note that the representation (J2]) 
makes sense for any q G R, although §T§ initially supposed that q G N. We shall also 
admit any value of v G R, although only v > — 1 corresponds to a real coupling, K G M. 

On a given lattice G, the Potts model can in general only be exactly solved along 
certain curves in the (q, v) plane. This is in sharp contrast with the Ising model (alias 
the q = 2 state Potts model), which can be solved at any temperature v [5]. Another 
major difference is that the Potts model has been solved only on a very few lattices, 
including the square [6] and triangular lattices [7J and certain decorated versions thereof 
[8]. On the contrary, the Ising model can be solved on essentially any lattice (provided 
that suitable boundary conditions, usually doubly periodic, are imposed). 

A first step in the solution of the Potts model on G is to determine the critical 
manifold — or critical frontier — , by which we understand the points in (q, v) space at 
which the model stands at a phase transition^ It is a remarkable fact that, at least 
in the solvable cases [6j QUI UJ, the loci of exact solvability coincide precisely with the 
critical manifold. It is equally remarkable that the critical manifolds on the square 
P, [TO] , triangular [7] and hexagonal (the dual of the triangular) lattices turn out to be 
given by simple algebraic curves: 

(f 2 — q) (v 2 + 4v + q) = , (square lattice) (3) 
v 3 + 3v 2 — q = , (triangular lattice) (4) 

v 3 — 3qv — q 2 = . (hexagonal lattice) (5) 

It is an outstanding question of lattice statistics to determine whether these features 
hold true more generally. 

The case of the Potts model on the kagome lattice has attracted particular attention. 
This is due not only to the practical applications of this model, but also to the fact that 

| The exact solution on the square lattice |6| shows that in the ferromagnetic regime (v > 0) the nature 
of the phase transition is first order for q > 4, and second order (continuous) for < q < 4. The belief 
that this is true on any lattice is supported by other exact solutions |7J and general field theoretical 
arguments. The nature of the transitions for q < and/or v < is obviously a more delicate question, 
which has been addressed only in the simplest cases [9] . 
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the kagome lattice is the "simplest" lattice on which an exact solution has not yet been 
found. More than thirty years ago Wu conjectured [11] that the critical manifold is 
given by the sixth-order algebraic curve 

v 6 + Qv 5 + 9v 4 - 2qv 3 - 12qv 2 - Qq 2 v - q 3 = . (kagome lattice) (6) 

The derivation of ([6]) relied on a number of exact equivalences and — crucially — a certain 
non-rigorous homogeneity assumption!^ 

The curve (Q has indeed a certain number of pleasing features: 

(i) For q — 2 it factorises as {v + l) 2 (w 4 + 4t> 2 — 8v — 8) =0. This situates the 
ferromagnetic phase transition at v c = a/3 + 2 a/3 — 1, in agreement with the exact 
solution [13\\i 

(ii) It passes through the origin (q, v) = (0, 0) with infinite slope, giving as expected a 
model of spanning trees (see [T5]). 

(iii) It behaves asymptotically as v 2 oc q for q ^> 1, ensuring as expected first-order 
coexistence between the phases dominated by the terms A = E and A = in (E}. 

Subsequent numerical results have however unambiguously established that is not 
an exact result: 

(i) For q — y 1 it places the bond percolation threshold at p c = = 0.524429 717 ■ ■ -, 
which is not correct but close enough that it took many years to be definitively 
ruled out [T7j. Recent numerics [2D] gives p c = 0.524404 978 (5). 

(ii) For q — 3 it gives v c = 1.876 269 208 ■ ■ -. This can be contrasted with the 
estimate |18| v c = 1.876 456 (40), obtained by analysing the low-temperature series 
for the magnetisation, susceptibility, zero-field partition function, and specific 
heat. Numerical diagonalisation of the transfer matrix provides the more accurate 
estimate [20] v c = 1.876 458(3). 

(iii) For q = 4 the prediction of ([6} is v c = 2.155 842 236 ■■ -. This cannot be 
discriminated by the less precise series result [H] v c = 2.1561(5), but is ruled 
out by the transfer matrix result [20] v c = 2.156 20 (5). 

(iv) Discrepancies have also been observed for large g > 4 [T9] . 

Below we shall present further conclusive numerical evidence against ([6]). Despite this 
numerical refutation of the conjecture ([6]), it is nevertheless remarkable; viewed as an 
approximation it is extraordinarily precise, the accuracy being of the order 10 -4 or 10 -5 
(see also [20J). 

§ For completeness we mention that Tsallis [T2j has proposed an alternative conjecture for the kagome- 
lattice critical manifold. This has however been shown definitely not to be the exact expression. 
Moreover it is less accurate than Wu's conjecture ([6]). 

|| The Ising model on the kagome and hexagonal lattices can be related by combining a star-triangle 
and a decoration- iteration transformation [14]. It follows that (1 + «kag) 2 = 2(1 + t>hcx) — 1 in our 
notations. Only one of the solutions to ([5]), Vhcx = 1 + v3j leads to real solutions for i>ka g , namely 
w kag = -1 ± \/3 + 2^3. 
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In view of this evidence it is tempting to try to improve on by systematically 
fitting the numerics to other low-order algebraic curves with reasonable (integer) 
coefficients. Such attempts [THl El] have however proved inconclusive. Another line of 
research is to extend the approximation (jSJ) to other (decorated) lattices of the kagome 
type [HI ED]. 

In this paper we present a general method for obtaining approximations to the 
critical manifold of the Potts models defined on any two-dimensional regular lattice 
G. The infinite lattice G is obtained by tiling two-dimensional space by a certain finite 
subgraph B C G that we shall refer to as the basis. The embedding of B in G determines 
exactly how G is obtained as a tiling by the motif B and will be defined precisely below. 
Corresponding to each embedded basis B we define a two-parameter graph polynomial 
Ps{q ) v) which is closely related to the Tutte polynomial |23j; the approximation to the 
critical manifold then reads simply 

P B (q,v) = 0. (7) 

The precision of the approximation can be systematically improved by increasing the 
size, in a sense to be made more precise in section of B. 

The definition of Ps{q, v ) conceals a number of remarkable features: 

(i) When not exact, the approximation (jZj) turns out to be very precise, even for the 
smallest possible choice of B. The analogue of flS} can thus be computed by hand 
rather easily for any regular lattice of interest. 

(ii) For some lattices G — or for the Ising model (q = 2) on any G — the polynomial 
Ps(q,v) factorises for any choice of B. The zero set of the recurrent factor then 
turns out to provide the exact critical manifold. In this sense, the computation 
of Ps{q,v) can be used to detect whether the Potts model on G might be exactly 
solvable. 

(iii) The determination of Pe(g,v) for larger bases is well suited for exact computer- 
assisted calculations. We shall pursue this point of view in section [3] below. 

Here we illustrate the general method and the above remarkable features for the 
cases when G is the square or the kagome lattice. 

The square lattice serves as a benchmark, since its critical manifold is known 
completely [6j (TOj E2l |9j . The result (j3J) is recovered by applying the general method 
to the simplest choice where the basis B consists of 4 edges. We present extensions to 
larger bases consisting of 8, 16 and 32 edges that all contain the exact result (E]) as a 
factor; the remaining factor gives information about the phase transitions at the Beraha 
numbers inside the antiferromagnetic regime [9]. 

In the case of the kagome lattice, the simplest case when B consists of 6 edges 
reproduces Wu's conjecture (E]). We then present extensions to larger bases consisting 
of 12, 24 and 36 edges. These extensions systematically improve the agreement of the 
predicted critical manifold with the existing numerical results. In addition we present 
improved numerical results for the critical manifold — both in the ferromagnetic (v > 0) 
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and antiferromagnetic (v < 0) regime — obtained by exact diagonalisation of the transfer 
matrix. Comparing this to the result from the 36-edge basis we find an accuracy of the 
order or 10~ 7 in the ferromagnetic regime (v > 0), and both qualitative and 

quantitative improvements on ([6]) in the antiferromagnetic regime (v < 0). 

The specialisation of these results to percolation (i.e., Pb(1, v)) has previously been 
reported by one of us [24]. Percolation polynomials for other lattices are reported in 

[251 m [27]. 

The paper is organised as follows. We begin the next section by reviewing the 
derivation of exact critical manifolds for three-terminal triangular-type lattices. We 
then show how these results may be used to define the graph polynomial Pb{o,v) on 
any regular lattice by the contraction-deletion algorithm and use this to rederive Wu's 
conjectures [TT] for the checkerboard and kagome critical manifolds. In section [3] we 
describe how this polynomial is calculated on larger bases using a computer program, 
and use the following sections to report polynomials computed for the square and kagome 
lattices on bases of up to 36 edges. In section [6] we support our conjecture that the 
polynomial approximations converge to the exact values by comparing our results with 
numerical calculations of critical curves. 

2. The graph polynomial 

The triangular-lattice Potts model with arbitrary interactions in up-pointing triangles is 
exactly solvable [29] . In particular, there exists a nice duality argument determining the 
critical manifold exactly. Since this is the starting point of our construction we begin 
by reviewing it. 

2.1. Triangular-lattice Potts model with interactions in up-pointing triangles 

Consider a Potts model on the triangular lattice in which all interactions occur inside 
up-pointing triangles, as shown in Fig. UK Let the three spins around an up-pointing 
triangle A be labelled in cyclic order as 01, cr 2 , cr 3 , starting from the spin in the lower- 
left corner (Fig. [lb). The most general form of the local Boltzmann weight on A123, 
compatible with the Sq permutational symmetry of the spins, takes the form 

W 12 3 = C + Ci<5 2 3 + c 2 5i 3 + C 3 5 12 + C 4 5i2 3 , (8) 

where we have introduced the short-hand notation 5ij = 8{a i: aj) and 8ijk = 
5(<Ji,(Tj)6(crj,(Tk) f° r the Kronecker delta functions. 

Each term in the sum can be represented graphically as follows: 

A A A A A 

Co Ci c 2 c 3 c 4 

Expand now the product over A 12 3 to form the partition function Z. This defines Z in 
terms of clusters on the hexagonal lattice G = (V,E) with local weights c p . Each cluster 
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(a) 



(b) 



Figure 1. a) Triangular lattice with arbitrary interactions in up-pointing triangles, 
b) Labelling of spins around an up-pointing triangle A123. 



comes with a weight q from the summation over the spins a, so that 

z = J2 * k[A) X\^) Np > ( 9 ) 

ACE p=0 

where N p is the number of up-triangles of type c p , and k(A) has the same meaning as 
in (ED. 

Alternatively, one can introduce loops that wrap around the boundaries and internal 
cycles of the clusters. The loops are represented graphically as follows: 






Go '0/2 <^3 'O4 



and a local weight a p is associated with each diagram. In down-pointing triangles, there 
is no interaction and each face is of the type ao (turned upside down) . One has therefore 
a model of loops on a triangular lattice — shifted vertically by 1/ \/3 lattice spacings with 
respect to the one on which the Potts spins are defined — and each edge is covered by a 
loop. 

Using the Euler relation the number of clusters k(A) and the number of loops 1(A) 
are related by 2k(A) = \V\ — \A\ + 1(A), where we recall that A is the set of blue edges 
on the hexagonal lattice G = (V,E). The partition function in the loop representation 
therefore reads 

z = ^'^n(a p y\ (10) 

loops p=0 

where n = ^Jq is the weight of a loop, and the local weights a p and c p are related by 

a =c , (11) 
naj = Q fori = 1,2, 3, (12) 
n 2 a 4 = c 4 . (13) 

Following [29] we now suppose that at criticality the model is invariant under a 
7r/3 rotation of the triangular lattice where the loops live. This implies that the critical 
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manifold is given by a 4 = a or [22] 

c 4 = gc . (14) 

Alternatively, the same relation can be found by a duality argument. Indeed, if we first 
transform the triangular lattice into a hexagonal lattice under duality, and next perform 
a decimation transformation of one half of the spins, we recover the original triangular 
lattice. Requiring that the weights be invariant under the combined transformation 
again leads to f lT4|) . 

2.2. Bases with three terminals 

As an example, we apply (|T4"]) to a triangular lattice with pure two-spin interactions, as 
in ([2]). We consider the case of arbitrary inhomogeneous two-spin couplings {^1,^2,^3} 
within A123, so that 

c = 1 , 

Ci = v i for i — 1, 2, 3 , 

c 4 = V1V2V3 + v\v 2 + v 2 v 3 + V3V1 . (15) 

The critical manifold is then given by ffT4"j) in the form ([7]), where we have associated 
with the three-edge basis B = A123 the graph polynomial 

P B (q, {vx,v 2 ,v 3 }) = v x v 2 v 3 + v x v 2 + v 2 v 3 + v 3 v x -q. (16) 

In the homogeneous case this indeed reduces to (jlj). 

We can also recover the special case of the square lattice by setting t> 3 = in ( TLB"!) . 
This formally corresponds to a two-edge basis (remove one edge from B), and in the 
homogeneous case one obtains the first factor in ([3]). 

It should be noted that the result f|T4l also applies when extra spins are present 
inside the up-pointing triangles^ One then simply sums out those spins to get an 
effective interaction of the form ([8]). As an example of this consider the hexagonal 
lattice with arbitrary inhomogeneous two-spin couplings {^1,^2,^3}. We have then 

c = vi + v 2 + v 3 + q , 
ci = v 2 ^3 , 
c 2 = v 3 v 1 , 
c 3 = ViV 2 , 

c 4 = v 1 v 2 v 3 . (17) 
The critical manifold again results from (fT4"j) in the form ([7j), with the graph polynomial 

P B (q, {vi,v 2 , v 3 }) = Wif 2 W3 - q(vi + v 2 + v 3 + q) . (18) 
In the homogeneous case we recover (JSJ). 

% This has been used extensively in [8]. One can also include explicit multi-spin interactions, but we 
shall not consider this possibility here. 
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Note that (ITS)) and (ITS!) are related by the duality transformation [2J v iV * = q, up 
to an unimportant global factor. More precisely 

Petri (q, {v*,vZ,v%}) = — P B he*(q,{vi,v 2 ,v 3 }) . (19) 

2.3. Square-lattice Potts model with checkerboard interactions 

Up to this point our determinations of the critical manifolds have been exact. We shall 
now see that crucial new ingredients appear when we consider larger bases. In particular 
we pay attention to the way the basis B is embedded in G so as to tile the entire lattice. 

Let us call a vertex of B a terminal if it will have to be identified with a vertex 
of one or more copies of B in the tiling of G. The remaining vertices of B are called 
internal. The models discussed in section 12.11 used a basis with 3 terminals, and in the 
embedding of B as up-pointing triangles each terminal was glued to a terminal in two 
other copies of B. This identification of each terminal with the other two is a convenient 
way of characterising the embedding. 

It should be clear from the examples ffl5|) and (ITT)) that the derivation of the 
weights Co and C4 are closely reminiscent of the computation of the partition function 
(T5J) itself. Indeed, these weights are just conditional probabilities that the terminals are 
connected in a certain way (not connected for Cq, and fully connected for C4) in the 
Fortuin-Kasteleyn expansion. 

The partition function ([2]) for the Potts mode|j] defined on any graph G = (V, E) 
can be computed by the contraction-deletion method (see [23]). Namely, let e G E be 
any edge in G. Denote by G/e the graph obtained from G by contracting e to a point 
and identifying the vertices at its end points (if they are different); and denote by G \ e 
the graph obtained from G by deleting e. The terms in the sum (J2J) can be grouped in 
two disjoint classes, according to whether e G A or e (ji A. For the terms with e G A 
(resp. e £ A), the number k(A) is unchanged upon contracting (resp. deleting) e. So the 
obvious generalisation of (j2J) to the case of arbitrary edge-dependent weights satisfies 

Z G (q, {v}) = v e Z G/e (q, {v}) + Z GV (q, {v}) . (20) 

Observing that Z factorises over the components of a disconnected graph, and that 
Z = q for an isolated vertex, the contraction-deletion formula (1201 allows one to compute 
Zc(q, {v}) recursively. 

We now define the graph polynomial Pb{q, {v}) corresponding to bases B with more 
than three terminals. It is computed recursively by applying the contraction-deletion 
formula (1201) to decrease the size of the basis. In the process we can rearrange B by 
identifying terminals from the embedding, and we shall avoid disconnecting B (modulo 
the embedding). The initial condition is to replace any 3-terminal basis obtained in the 
recursion by its corresponding critical manifold, which has been obtained in section 12.21 

+ In the mathematics literature, the Potts model partition function is known as the multivariate Tutte 
polynomial |23j . 
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(a) 



(b) 



Figure 2. a) The square lattice with the checkerboard couplings shown in b). 



Our key assumption — that we shall test in the following — is that P B (q,{v}) = 
provides an approximation to the critical manifold that becomes more and more precise 
(when it is not exact) as the size of B increases. 

Let us illustrate the above definition on a square-lattice Potts model with 
interactions in alternating faces forming a checkerboard pattern; see Fig. [2^. The 
corresponding basis has four terminals, and the embedding amounts to gluing 
diametrically opposite terminals. We are interested in the fully inhomogeneous case 
with couplings {vi, v 2 , v 3 , 1*4}, as shown in Fig. [2b. Performing the deletion-contraction 
of the edge with weight V4 we obtain 

P B (q, {vi,v 2 ,v 3 ,V4}) = v 4 P B tn(q, {vx,v 2 ,v 3 }) + P B ^(q, {fi,u 2 ,f 3 }) ■ (21) 

In the first term the two terminals have been identified, so we arrive at the known case 
of the triangular lattice. In the second term we have made use of the embedding to flip 
one of the edges. In doing that, one of the terminals is turned into an internal vertex, 
and we recognise the 3-terminal basis of the hexagonal lattice. Substituting now ( |T6l) 
and dUD into (J2TJ) we find 

P B (q, {V 1, V 2 , V 3 , Vi}) = V X V 2 V 3 V 4 + {v 2 V 3 V 4 + Vl?^4 + V X V 2 V 4 + V1V2V3) 

- q(vi + v 2 + v 3 + v 4 ) - q 2 . (22) 

Note that this expression has an S4 symmetry under the permutation of any two 
couplings. 

The expression P B (q, {v}) = with ( |22l) was first derived by Wu [H] using a 
different method (and a homogeneity assumption). He initially conjectured that it was 
the exact critical manifold. However, Enting (30] has observed that the putative S4 
symmetry is broken — albeit only at high order — in the series expansions for various 
physical quantities of the checkerboard-lattice Potts model, and he concluded that (|22|) 
cannot be correct [31] . In section 6.2 of [32] Maillard seems nevertheless to suggests that 
( |22|) may be correct. Unfortunately we do not know of any numerical investigation of the 
critical manifold for the completely inhomogenious checkerboard Potts model that would 
allow us to assess whether ( 1221) is correct, and if not, how accurate an approximation it 
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is. But in section @] below we have computed the polynomials Pb{q, {v}) corresponding 
to larger bases for the square lattice, with four different couplings arranged in a 
checkerboard pattern. We observe that these polynomials invariably factorise, shedding 
the small factor (f22|) . This is a strong indication that (122)) is indeed the exact critical 
manifold of the checkerboard-lattice Potts model. 

Something interesting happens if we specialise to the usual square lattice with 
horizontal couplings V\ = v 3 and vertical couplings v 2 = ^4- Then (1221) factorises as 

P B (q, {vi,v 2 ,v 1 ,v 2 }) = (viv 2 - q){viv 2 + 2v x + 2v 2 + q) , (23) 

and in the homogeneous case we recover Q. The zero set of (1231) is known to be 
the exact critical manifold. Indeed, the first factor corresponds to the ferromagnetic 
transition curve [6] and the second one to the antiferromagnetic transition [10] (see 
also [22] [9]). It is particularly remarkable that the antiferromagnetic transition curve 
is obtained exactly by this approach, since — to our knowledge — it cannot be obtained 
by a simple duality argument. Indeed, [10] relies on the solution of the Yang-Baxter 
equation and a so-called Z-invariance. 

2.4- Wu's conjecture recovered 

We now turn to the kagome lattice G. The simplest choice of basis B e is the 4-terminal 
arrangement of six edges shown in Fig. [3k which produces G by the checkerboard 
embedding discussed in section 12.31 Using the contraction-deletion identity (120]) on 
the t>4 edge we obtain for the corresponding graph polynomial K({v}) = PB 6 {q, {v})'- 

K(v 1 ,v 2 ,v 3 ,v^,v 5 ,vq) = ViB(v b + v 6 + v 5 v 6 ,v 1 ,v 2 ,v 3 ) 

+ A(v 1 ,v 2 ,v 3 ,v 5 ,v 6 ) . (24) 

In the first term, the w 4 edge has been contracted, so that the v 5 and the v e edges are 
in parallel. Using parallel reduction [23J these can then be replaced by a single edge 
with coupling v 5 + v 6 + v^vq. The result is a 3-terminal basis of the martini-B lattice 
[33| 134"] (see Fig. |3b) whose graph polynomial we have denoted B(vi,v 2 , t> 3 , t> 4 ). This is 
then replaced by the critical manifold ( [141 which reads explicitly [35] 

B(Vl, V 2 , V 3 , Vi) = Wi(w 2 ^3 + V 2 V 4 + + v 2 v 3 v 4 ) 

- q{q + v 1 + v 3 + v A ) . (25) 

In the second term in ( 124]) the t> 4 edge has been deleted. Flipping the v 5 and v q edges 
using the embedding, the result is a 3-terminal basis of the martini-A lattice (see Fig. |3b) 
whose graph polynomial we have denoted A(v\, v 2 , v 3 , v±, V5). Using again (114]) it is 
replaced by the critical manifold [35] 

A(vx, V 2 , V 3 , Vi, V 5 ) = V4V5(VlV 2 + V 2 V 3 + V 3 V! + V!V 2 V 3 ) 

- q(q 2 + q(vi + ^2 + ^3 + ^4 + v 5 ) + ViV 2 v 3 

+ (v 2 + v 4 )(v 3 + v 5 ) + vi(v 2 + ^3 + ^4 + ^5))- (26) 
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Inserting ( 125|) and (126|) into (|24|) produces the desired graph polynomial. The 
corresponding approximation for the critical manifold is K({v}) = by (|7|). The case 
of a kagome lattice with different couplings for each of the three principal directions 
corresponds to the choice v± — Vi, v$ — t>2 and vq = V3. This coincides with Wu's most 
general conjecture for the critical manifold (see Eq. (11) in [TTJ) and reduces to in 
the homogeneous case. 

We have already discussed in the introduction that (JSJ) is a very precise 
approximation to the critical manifold, but not an exact result. Further numerical 
evidence for this statement will be given in section |6j 

3. General bases 

As remarked above, the general contraction-deletion formula for Pb{q, {v}) is identical 
to that for the partition function (see f !2Up ). namely 

P B (q, {v}) = v e P B/e (q, {v}) + P B \ e (q, {v}) . (27) 

In the simple kagome example given above, both contraction and deletion yielded lattices 
for which the exact solution is known. For larger bases, this will not always be the case 
and then the algorithm must be applied to whichever (or both) of B/e and B \ e is 
unsolved, and so on recursively until solved lattices appear. This results in a sort of 
binary tree in which each node is a lattice and branches terminate on known graphs. 
However, the final polynomial is independent of the edge chosen at each step, and is 
therefore a unique property of the graph B and the terminal identifications that define its 
embedding in G. Of course, the complexity of the problem increases exponentially with 
the number of edges in the basis, and thus the algorithm must be handled by computer 
for large bases. The operation of the program that implements the contraction-deletion 
algorithm can be sketched as follows: 
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(i) As input, it takes an array of internal vertices, terminals and edges, along with 
the specification of their connectivities and identifications between terminals that 
determine how the basis is tiled to form the lattice. 

(ii) An edge, e, is then chosen, respecting a limited set of criteria (e.g., not disconnecting 
the basis), for contraction and deletion to give the new bases G/e and G\e. 

(iii) Often, these new graphs have edges doubled in series, or, as for the martini-B lattice 
in $M§, in parallel. These may be simplified by replacing them by a single edge 
with an effective interaction. Other complications are also possible, such as isolated 
or dead-end vertices, and these are removed from the problem. 

(iv) After simplifying, G/e and G\e are checked against a small set of known lattices. 
These are the square, hexagonal and triangular lattices, along with the one- 
dimensional chain and certain situations that are known to give critical manifolds 
that are identically zero (see [2H])- We restrict ourselves to these lattices because of 
the relative ease of identification (e.g., all edges are equivalent and we are thus not 
required to worry about edge-matching). If a graph is recognised, it is replaced by 
the corresponding graph polynomial and the process ends. If not, step 2 is called 
recursively with the unknown graph as input. 

(v) The final result is a tree of lattices with branches that terminate on known 
lattices. The output of the program is a list of functions giving the contraction- 
deletion formula for each lattice. This list is evaluated in Mathematica and the 
homogeneous polynomial can then be found. 

This program was used to compute percolation (q = 1) polynomials in [21] and [28J. 
Only minor adaptations were required to find the polynomials for general q, so we discuss 
only those changes here and refer the interested reader to [28J where the implementation 
is described in detail. 

3.1. Simplification 

The only differences between the percolation and Potts implementations are in the 
simplification phase of the algorithm, most importantly in the way interactions are re- 
labelled upon the removal of doubled bonds. The handling of dead-end vertices also 
requires a minor modification. 

3.1.1. Edges doubled in parallel. Contracting the V\ edge in the martini-A lattice 
of Fig. [3b leaves t> 2 and v 3 doubled in parallel. These must be replaced by a single 
edge. In the percolation computation, the probabilities P2 and p 3 are combined into 
1 — (1 — j9 2 )(l — ps) = P2+P3 —P2P3, i-e., the probability that at least one bond is open. 
The Potts expression for the effective interaction, v p , is similar [23J: 



v p = v 2 + v 3 + v 2 v 3 . 
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3.1.2. Edges doubled in series. For the martini-A lattice in Fig. [3b, deleting the v\ 
edge leaves edges doubled in series, for example v 3 and v 5 . These are replaced by the 
effective interaction, v s , given by |23j : 

(29) 



v 3 + v 5 + q 

The denominator in (j29p presents a slight complication. Consider the deletion of v\ in 
Fig. [3b, giving the square lattice, which has the inhomogeneous manifold S(vi,v 2 ) = 
V\V2 — q = 0, but with edges doubled in series. The contraction-deletion formula 
( |27l) works because the completely inhomogeneous polynomial is at most first-order 
in any interaction; we may have terms like V1V2V3 but not vfv 2 or, more to the point, 
V\/{v2 + V3 + q). As such, we cannot use 

S( V3V5 , V2V * ) (30) 



v 3 + v 5 + q V2 + V4 + q, 

for the manifold of B \ e. However, because the critical manifold is found by setting 
S(v\,V2) = 0, we are free to multiply away these denominators, and we must do this 
before inserting an expression into the contraction-deletion formula. Thus, we have 

Mpx, v 2 , v 3 , u 4 , v 5 ) = v 1 H(v 2 + v 3 + V 2 V 3 , V 4 , V 5) 

+ (.3 + v 5 + q)(v2 + v 4 + q)S ( V J V \ , V f\ ) , (31) 

from which we recover (1261) . 

3.1.3. Dead-end vertices. An internal vertex connected to only one edge is a dead end. 
In percolation, such an edge contributes nothing to the connectivities between terminals, 
and therefore nothing to the graph polynomial. Its probability, p, disappears from the 
problem because it only appears in an overall factor of [p+ (1 — p)] = 1 multiplying the 
probability of every event. So the edge is simply removed from the problem along with 
the vertex. In Potts language, the weight of every event is multiplied by (v + q) and 
thus, although we still remove the edge and vertex, the manifold must be multiplied by 
this factor upon the removal. We give a simple example. In Fig. Hk, we have a triangular 
lattice with two edges doubled in series. Adhering to the rules set out in the previous 
section, the graph polynomial for this is given by 

(v3 + v A + q)T(v 1 ,V2, V3V \ V (32) 



V3 + v 4 + q, 

Contraction-deletion should also produce this answer. Contracting t> 4 just gives the 
simple triangular lattice, T(%, v 2 , v 3 ), but deleting this edge results in the square lattice 
with a dead-end vertex and edge, v 3 . We remove these from the problem but retain the 
factor of (q + v 3 ) to give (q + v 3 )S(v 1, t> 2 ). Injecting these into the contraction-deletion 
formula, 

v 4 T(v 1 ,v 2 ,v 3 ) + (q + v 3 )S(v 1 ,v 2 ), (33) 
we do indeed recover (1521 . which can be seen by expanding both expressions. 
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Figure 4. a) Triangular lattice with a series-doubled edge; b) deletion of the V4 edge 
leaves V3 and its right vertex as a dead end. 
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Figure 5. a) 8-edge; b) 16-edge; and c) 32-edge bases for the square lattice. Terminals 
with matching shapes are identified. 



Aside from these changes, the main operation of the program, including edge 
selection and lattice identification, are exactly as described in [28J. We also have the 
same upper limit of practical feasibility, namely using basis with at most 36 edges. 

4. Square-lattice Potts model using larger bases 

Summarising, we have seen that the definition of Pb(q,v) made in section [2731 leads to 
known approximations of the critical manifold when the lattice G is obtained from a 4- 
terminal basis. In the remainder of the paper we shall corroborate our key assumption 
that by using larger bases the accuracy of the approximations can be improved, and 
more subtle features of the critical manifolds can be uncovered. 

To this end we study first the benchmark example of the square lattice in this 
section. The following section will be devoted to the kagome lattice. 

We first extend the 4-edge checkerboard to the 8-edge basis in Fig. [5k- The 
polynomial for this case is given by 

(q - v 2 ) (q + Av + v 2 ) (q 2 + Aqv + v 4 + 4w 3 + 8v 2 ) , (34) 

and we see the exact solutions appearing as factors. The polynomial for the 16-edge 
basis depicted in Fig. [5b is given in the Appendix in equation flAlj) . We also found the 
polynomial for the 32-edge basis in Fig. |5b, which is given in equation ( 1A2I) . 
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Figure 6. Manifolds Psiq, v) for various choices of bases B for the square lattice (see 
text). The exact critical manifolds ([3]) obtained from the bases with 2 and 4 edges are 
contained in the manifolds corresponding to larger bases. 



We notice for the polynomials (152]) . (IXIji and (1X2]) that 

P B (g,i;) = -t;l £ lg-l £; l/ 2 P B (g,^) (35) 

with f * = g/f , where \E\ is the number of edges in the basis. This means that the graph 
polynomials respect the self-duality of the square lattice. In particular, the solutions to 
Pb{q,v) = are either self-dual curves, or pairs of mutually dual curves. 

It appears that the exact solvability [6j [10] of the square-lattice Potts model 
manifests itself in the presence of the factors fl3]) in all the graph polynomials constructed 
from larger bases, irrespective of their embedding. More generally, we have found that 
for four different coupling constants arranged in a checkerboard pattern (see Fig. |2]), the 
polynomials always contain the factor ( 122]) . It remains to assess whether the remaining 
(in general non-factorisable) factor in Pg has any physical relevance. 

To investigate this issue, we have plotted the manifolds Pb{q, v ) = for the different 
bases in Fig. [6j It is seen that the remaining factor produces dual pairs of curves inside 
the dual pair of antiferromagnetic transition curves v 2 + 4t> + q = 0, in the form of 
"bubbles" emanating from the point (q,v) = (4,-2). Each bubble intersects the self- 
dual transition curve v = —^/q in two points, of which one is exactly go — 4 and the 
other we denote q c . The bubbles are almost vertical in the vicinity of that latter point. 
For the 32-edge basis we find two bubbles with respectively 

q cA = 2.000 059 024- •• , 

g Cj2 = 3.088 542 134 • • • . (36) 
These values are conspicuously close to 2 and 3. 
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It was shown in [22J E] that the Potts model possesses singularities inside the so- 
called Berker-Kadanoff phase, which is the region bordered by the antiferromagnetic 
transition curves v 2 + 4t> + q = 0. These singularities occur at the Beraha numbers 
q = Sfc with 

B k = (2cos(tt/A;)) 2 , with k = 2, 3, 4, . . . (37) 

and are independent of v. For a given integer k the singularity may or may not concern 
the dominant term in the partition function, and only the former case leads to a phase 
transition. The independence of v implies that the incidence on the phase diagram is 
the formation of a vertical ray in the (q, v) plane. 

The issue of dominance depends on the choice of boundary conditions. In [36] 
it was argued from results of conformal field theory — and checked numerically — that 
with cyclic boundary conditions (free in one lattice direction and periodic in the other) 
vertical rays are formed for k e 2N. The corresponding result for toroidal boundary 
conditions [37] (periodic in both lattice direction) is that vertical rays occur only for 
k — 4 and k — 6. Because of the identification of opposite terminals in the embeddings 
of the basis B that we have used, it is the toroidal boundary conditions that are relevant 
in the present case. 

In (|36|) we have indeed q c ^ ~ £> 4 = 2 and g Cj2 ~ B 6 = 3, in agreement with [37]. We 
therefore conjecture that in the limit of an infinitely large basis, PB{q,v) will contain 
(EJ) and a remaining factor whose zeros produce a couple of vertical rays at q = 2 and 
q = 3, extending between the two branches of the antiferromagnetic transition curve 
v 2 +4v + q = 0. The two bubbles obtained from the 32-edge basis provide strong support 
of this conjecture; we believe that the non-vertical parts of those bubbles will coincide 
with parts of the antiferromagnetic transition curve in the limit of infinite B. 

To produce an almost vertical ray from Ps(q,v) = clearly requires the algebraic 
expression Pe(g, v) to be of large degree in q and v. Obviously the degree in q equals the 
number of vertices in the basis B (counted up to identification of terminals through the 
embedding), while the degree in v equals the number of edges in B. The above conjecture 
is compatible with the increase of degree as the basis becomes larger: infinite-degree 
expressions are needed to produce truly vertical rays in the limit of an infinite basis. 

5. Kagome-lattice Potts model using larger bases 

We are now ready to fulfill the main objective of this paper, which is to improve on Wu's 
approximation ([6]) for the critical manifold of the Potts model on the kagome lattice. 
Starting with the 12-edge basis in Fig. [7_H, we find the polynomial 

(q 3 + 6q 2 v + Ylqv 2 + 2qv 3 - 9v 4 - 6v 5 - v 6 ) x 

{q 3 + 6q 2 v + 16qv 2 + 2Av 3 + 2qv 3 + 17v 4 + Qv 5 + v 6 ) = . (38) 

The first term in brackets is once again Wu's formula ([6]), so, similar to the situation 
in percolation |24j, we get no better estimate of the ferromagnetic critical point by 
extending from 6 to 12 edges. However, the second term does provide us with additional 




Figure 7. a) 12-edge; b) 24-edge; and c), d) 36-edge bases for the kagome lattice. 



information in the antiferromagnetic region, which is plotted in red in Fig. [HJ We also 
found the polynomials for the 24-edge basis of Fig. [7b, and the two different 36-edge 
bases of Fig. [Tfc and [7U. These polynomials are given in the Appendix (see eqs. flA3j) . 
( ED , and flMp ). 

The manifolds Ps{q,v) = for the different bases are plotted in Fig. [HJ 

In the ferromagnetic region (v > 0) the difference between the various curves is 

not visible on the scale of the figure. The approximations for the critical coupling in 

percolation read 

{ 1.102 738 621- ■■ (6-edge basis) 
1.102 636 956- •• (24-edge basis) (39) 
1.102 632 538 •• • (36-edge basis) 

to be compared with the most precise numerical estimate [20] v c (q = 1) = 
1.102 629 24(2). The corresponding percolation threshold is p c = j^j-- For the Ising 

model (q = 2) all curves pass through the exact critical point [13] v c = a/3 + 2\/3 — 1 = 
1.542 459 756 • • •. When q = 3 we have 

{ 1.876 269 208 •• • (6-edge basis) 
1.876 439 754- •• (24-edge basis) (40) 
1.876 447 147 •■ ■ (36-edge basis) 

to be compared with the numerical estimate |20j v c (q = 3) = 1.876 458 (3). Finally, for 
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Figure 8. Manifolds Pe(g,w) = for various choices of bases B for the kagome 
lattice (see text). The manifold corresponding to the basis with 6 edges is contained 
as a factor in the 12-edge manifold. The letters in the figure serve as a guide to the 
discussion in the main text. 



q = 4 we find 

{ 2.155 842 236 •• • (6-edge basis) 

2.156 207452- •• (24-edge basis) (41) 

2.156 223 187 ■■ ■ (36-edge basis) 

whereas the numerics |20j gives v c (q = 4) = 2.156 20 (5). 

We defer a more detailed discussion of the ferromagnetic region to section El 
where we shall also present our own numerical results. Instead we now discuss some 
common features and differences of the curves and formulate some conjectures about 
the thermodynamical limit. 

There is a branch that extends to infinite q in both the ferromagnetic and 
antiferromagnetic (v < 0) regions. Both of these have the correct asymptotic behaviour 
v oc y/q for all bases. Indeed, on a lattice of coordination number z, one needs v oc q z l 2 
in order to ensure first-order phase coexistence between the two dominant terms (with 
A = E and A = respectively) in the expansion t$Z§. 

With all bases, the ferromagnetic curve runs into the origin O : (q,v) = (0,0) with 
infinite slope, as is required [12] to produce a model of spanning trees. The continuation 
into the antiferromagnetic region bends around at A : (q,v) ~ (0.96, —1.58), and then 
goes down to the point B : (q, v) = (0, —3). Since all curves pass through (q, v) = (0, —3) 
exactly, we conjecture that this is a correct feature of the thermodynamical limit. If 
correct, this would imply that on the dual (diced) lattice, the problem of spanning 
forests pj)] has a critical point with a weight per tree w c = —3. 
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After bending around in B the curves continue in the region v < —3 to a point 
C : (g, v) m (4.35, —3.03) where they bend around again. It is interesting that C 
definitely has q > 4, a feature that was not present in Wu's conjecture (alias the 6-edge 
curve). 

The curves show strong finite-size effects in the domain —3 < v < (excluding the 
branch OAB which has been discussed above). In particular, it might be that the point 
C will eventually move to larger q, or even infinity. 

For all bases the curves pass through D : (q, v) = (2, —1) exactly, and we conjecture 
that this is a correct feature of the thermodynamical limit. For larger bases the branches 
BD and DC have a tendency to close up for q w 2 so as to form an almost vertical ray. 
We conjecture that in the thermodynamical limit there will be an exactly vertical ray 
extending between E : (q,v) = (2,t>i) and F : (q,v) = (2,t> 2 ), for some v\ < t>20 

It in then natural to assume that BEC and OFoo will be smooth transition curves 
in the thermodynamical limit. The 36-edge basis shows that the space in between 
these two curves is likely to contain a further vertical ray. Since the properties of 
the Berker-Kadanoff phase can be argued to be universal [39] we conjecture that the 
other ray will again be located at q = 3 [37], cf. our discussion of the square lattice 
above. We further conjecture that the curves OAB, BEC and OFoo will surround 
the Berker-Kadanoff phase. To delimit this phase to the right, one more curve is 
needed. The 36-edge basis provides a convincing finite-size estimate for this latter 
curve (starting at (q, v) ~ (3.06,-1.92), bending around at (q,v) ~ (3.32,-1.39) and 
ending at (q,v) « (3.22,-1.08)). 

We believe that in the thermodynamical limit the curve Foo will pass through the 
point G : (q,v) = (3, —1) exactly. Namely, one can show [40] that the three-state zero- 
temperature antiferromagnet on the kagome lattice is equivalent to the corresponding 
four-state model on the triangular lattice. The latter is known to be critical with 
central charge c = 2 (see [40] and references therein). Our finite bases locate the 
antiferromagnetic transition in the q = 3 model at 



and it seems likely that this might tend to v^ F (q = 3) = — 1 in the thermodynamical 
limit. This is presumably the point in the phase diagram where the finite-size effects 
are the most important. 

We should stress here that to get a reliable picture of the phase diagram it is not 
sufficient to study just one basis B, however large. Indeed it is obvious from Fig. [8] that 
different bases reveal different parts of the critical manifold. For instance, the 6-edge 

* Note that the existence of this vertical ray is not in contradiction with the argument of Huse and 
Rutcnberg [38] (using results of [14]) that the kagome-lattice Ising model is disordered for any v > — 1. 
Indeed, the difference between sitting at q = 2 exactly, and taking the limit q — Y 2 in the (q, v)- 
plane phase diagram, is the quintessence of the special physics at the Beraha numbers within the 
antiferromagnetic region [22| [9] . 




0.921 400 117 •• • (6-edge basis) 
0.973 665 377 ■• • (24-edge basis) 



(42) 
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basis misses completely the curve OFEB. And even though the 36-edge basis gives in 
many respects the most precise approximation to the true critical manifold, it misses a 
part of the curve Foo which is covered by the smaller 24-edge basis. 

Quite obviously it would require more work to give numerical support for the phase 
diagram discussed above and determine the relevant critical properties. In section [6] we 
limit ourselves to a detailed investigation of the ferromagnetic transition curve and of 
the antiferromagnetic branch BC. 

Note finally that the slope 1/w with which the curve OD goes into the origin 
determines the critical point in a model of spanning forests [15J on the kagome lattice. 
The critical weight per component tree can be determined as 

{-3.364 655 607 ■• • (12-edge basis) 
-3.553 344 713- •• (24-edge basis) (43) 
-3.578 781346- •• (36-edge basis) 

We can also determine the critical point of the flow polynomial which is 

related to the Potts- model partition function fl2]) by setting v = — q [23]. We find 

{ 3.324 717 957 ■• • (6-edge basis) 
3.400 923 464- •• (24-edge basis) (44) 
3.405 701 476 •■ ■ (36-edge basis) 

By duality this can also be interpreted as the critical point of the chromatic polynomial 
on the diced lattice. Our results (144"]) agree with the crude extrapolation q c m 3.4 given 
in [49]. 

6. Numerical results 

We have performed extensive numerical simulations in order to accurately locate (parts 
of) the critical manifold for the kagome-lattice Potts model. The objective is to assess 
whether the graph polynomials found in section |5] improve on Wu's conjecture (EJ) — and 
if so, by how much — when the size of the basis increases. 

For a critical system, the finite-size scaling of the free energy per unit area f(L) on 
a cylinder of circumference L reads j4TJ [42] 

TVf* 

/(L) = /(oo)- — + o(L- 2 ), (45) 

where /(oo) is the free energy in the thermodynamical limit and c is the central charge. 
This result also provides an accurate means of locating the critical point, since the 
numerically determined c will exhibit a local extremum at criticality [43J. The precision 
can be improved [44j by adding a non- universal term A/L 4 to (145!) . 

6.1. Transfer matrix construction 

The quantity f(L) can be obtained from the largest eigenvalue of the transfer matrix. 
There are two natural transfer directions for the kagome lattice: parallel or perpendicular 
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to one third of the lattice edges. In Fig. [7J the parallel and perpendicular transfer 
directions are respectively horizontal and vertical. Let T^v be the transfer matrix written 
in a basis of states that keeps track of N spins in a "time slice" — a layer perpendicular to 
the transfer direction. Using the basis provided by the Fortuin-Kasteleyn representation 
02]) then determines the number of basis states (i.e., the dimension of T/v) as the Catalan 
number jl5] 

c " = 7vTl(T)-]wi 1 + 0(1/7V) l' < 46 » 

The physical width L is N times the height (resp. N times the side length) of an 
elementary triangle for the parallel (resp. perpendicular) transfer direction. It is 
therefore most efficient to choose the perpendicular direction (i.e., such that the time 
runs vertically in Fig. [7]). Setting the side length of an elementary triangle to unity, we 
can henceforth identify N = L. 

We suppose henceforth L even, since otherwise the kagome lattice would not be 
compatible with periodic boundary conditions across the time slice. Let I, J^+i and 
Dj be the identity, join and detach operators (see [16]) satisfying the Temperley-Lieb 
algebra [UJ with weight q per connected component. The composite operators 

Hi = I + t>J M+1 , 

V, = v\ + D, (47) 

then add respectively a horizontal ("space-like") edge between the vertices i and i + 1, 
and a vertical ("time-like") edge a vertex i. The edge weight v is that of ([2]). From these 
we can define an operator that builds the basic bow tie pattern of the kagome lattice: 

B; = HjVjHjDj+iHjVjHj . (48) 

The transfer matrix then reads (by simple inspection of Fig. [7J 



Tl = I T B* TT B 2-i • (49) 




Since this propagates the time slice by four times the height of an elementary triangle, 
the corresponding free energy per unit area is 

/(L) = -^logA max , (50) 

where A max is the eigenvalue of Tl with maximal norm. 

We have diagonalised Tl for L = 2,4,6,.. .,16 and computed finite-size estimates 
c(L) of the central charge c from ( I4"5j) with the A/L A term included, by using three 
successive sizes (L — 4, L — 2, L). For each fixed q = ^,^,...,|[] we have varied v 
until a local extremum was found. This extremum provides an estimate of the critical 
coupling J^| 

jj For q = 1 there is no finite-size dependence in /(£), so the method does not work in that case. 
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Figure 9. Estimates of the critical curve in the ferromagnetic regime, relative to the 
approximation provided by the 24-edge basis. The solid lines show the zeros of Pb(q, v) 
for various bases B. The points are numerical transfer matrix results for system size L 
(see text). Extrapolations for L — > oo are shown for q = 3 and q = 4. Other numerical 
determinations are given for q = 1, q = 3, and q = 4 (and labelled "other work"). 



6.2. Ferromagnetic critical curve 

We first present the results in the ferromagnetic region v > 0. It is evident from Fig. [8] 
that the various approximations to the critical manifold are extremely close, so in order 
to appreciate the differences we focus on the quantity v — 1>24, where t>24 is the relevant 
root of Ps(q,v) using the 24-edge basis. 

Fig.[9]shows our numerical results (up to size L = 16) along with the approximations 
obtained from the graph polynomials of section [5] using bases with 6, 24 and 36 edges. 
For q = 3 and q = 4 we further show extrapolations of the numerical results: 

1.876 459(2) (q = 3) 

2.156 252(2) (q = 4) 1 ' 

These were obtained by fitting the residual dependence in 1/L to first and second-degree 
polynomials, and estimating the error bars by comparison and successive elimination of 
the data points with small L. We also show in the figure the previous numerical results 
for q — 1 (percolation [TBI EOl); 9 = 3 (series expansions [IS]), and q = 4 (transfer 
matrix diagonalisations [20j) lf fl 

It is clear that at this level of precision the Wu conjecture (JBD — i.e., the 6-edge basis 
result — is definitively invalidated. More precisely, it is ruled out by our results (T5TT) with 

ff We do not show the series result [TH] for q = 4, since the corresponding error bar is several times the 
vertical extent of Fig. [9] The q — 3 result of [20] is also not shown, because it is almost identical to 
our own extrapolation. 
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a confidence level of 95 (resp. 205) standard deviations for q = 3 (resp. q = 4). 

More importantly, it is obvious from Fig. [9] that the improved approximations with 
24 and 36-edge bases move systematically towards the numerical results. The results 
from the 36-edge base are at a tiny distance from the numerics, of the order 10 -6 or 
10~ 7 . It is nevertheless clear that even those results only constitute an approximation. 

Note also that the numerics provides (within machine precision) the exact v c for 
the Ising model (q — 2). 

We finally remark that in a recent paper Baek et al. [IE] located the critical points 
for q = 3 and q = 4 using phenomenological scaling of the internal energy crossings. 
These we obtained from transfer matrices (with the "parallel transfer direction" in 
the terminology of section I6.ip . The authors of [48] claimed that these results are 
independent of system size, whence the results for the smallest possible systems (N = 2 
and N = 4 in the notation of section IBTTl) : 

v c = 1.876 313463 895- •• , (9 = 3) 

v c = 2.156174166 284- •• , (q = 4) (52) 

would be exact. This claim is however not true. Internal energy crossings indeed 
exactly determine the critical point for the square-lattice Potts model (by an easy duality 
argument), and for the Ising model on more general lattices including kagome. However, 
for the q ^ 2 Potts model on the kagome lattice, the energy crossings do exhibit finite- 
size corrections. Note also that the results ( 152]) are incompatible with our numerical 
results ( l5il . 



6.3. Antiferromagnetic critical curve 

We have similarly investigated numerically the antiferromagnetic critical curve in the 
region v < —3 (referred to as BC in section [5]). As in the ferromagnetic case we show 
the differences v — v 2 4 (see Fig. ITU]) . 

It is evident from the scale of the vertical axis that the finite-size effects are much 
larger than in the ferromagnetic region. In particular, larger bases would be needed to 
convincingly approximate the numerical results, especially in the region q ~ 4. 

One exception is again the Ising model (q = 2) where finite-size effects are 
completely absent. Indeed, all curves pass through the exact value [13] v c (q = 2) = 
-V3T2V3- 1 = -3.542 459 756 

Extrapolations of the finite-size numerics — obtained using the method described 
above — read: 

'-3.486(2) (g = 3) 

-3.361 (14) (q = 4) K ' 

and are shown in Fig. [101 The result for q = 3 is in agreement with the Monte Carlo 
determination [49] of the critical coupling of the 3-state Potts model on the diced lattice, 
^diced = _o.860 599(4). The corresponding dual value v c = 3/vf ccd = -3.485 94(2) is 
in fact more precise than (1531) . 
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Figure 10. Estimates of the critical curve BC in the antifcrromagnetic regime, relative 
to the approximation provided by the 24-edge basis. The symbols and colours have 
the same meaning as in Fig. [5] 



The Monte Carlo study (H] also shows that the 4-state Potts model on the diced 
lattice is disordered for all — 1 < -y dlccd < 0. The corresponding dual statement is that 
for q — 4, there is no phase transition for — oo < v < —4. This again agrees with the 
fact that our critical manifolds avoid this interval. 

6.4- Critical properties 

The numerical results for the central charge along the ferromagnetic transition curve 
are shown in Fig. [HI They agree with the universality class of the usual ferromagnetic 
Potts model: 

- 1 -M^T)' < 54 > 

where q = Bk is parameterised by k G [2, oo) through fl37j) . In fact, the deviations from 
the exact result are not discernible on the scale of the figure. 

The corresponding results along the antiferromagnetic curve BC are given in Fig. [121 
Remarkably they again agree with (1541) . Note that there are now rather strong finite-size 
effects in the region q ~ 4. 

The result for k = 6 — namely that the q = 3 model at the antiferromagnetic 
transition ( [53]) is in the universality class of the three-state ferromagnetic model — was 
previously reported in 
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7. Discussion 

In this article we have defined a graph polynomial Pb{q,v) that characterises a finite 
basis graph B and the way it is embedded in order to tile an infinite regular lattice 
G. Being closely related to the Tutte polynomial, Pb((7, v) can be evaluated recursively 
from the contraction-deletion formula ( 120|) . 

The zero sets ([7]) were shown to reproduce well-known approximations for the 
critical manifolds for the Potts model on the checkerboard and kagome lattices [H] 
in the case of the smallest possible bases. We have shown that such approximations 
can be systematically ameliorated by increasing the size of the basis. In particular, the 
36-edge basis for the kagome lattice provides an approximation of the critical manifold 
whose deviation from our high-precision numerical study is of the order 10~ 6 or 10 -7 
in the ferromagnetic regime (v > 0). In the antiferromagnetic regime (v < 0) our 
approximations reveal qualitatively new structures, such as singularities at the Beraha 
numbers ( 1371) and the extent of the Berker-Kadanoff phase [221 |9] . 

The fact that the resulting 32nd-order algebraic curve still does not reproduce 
perfectly the numerical data leads us to conjecture that the true critical manifold for 
the kagome lattice is non-algebraic. It is tempting to speculate that this behaviour may 
be responsible for the fact that the Potts model on the kagome lattice has not yet been 
solved (except for the q = 2 Ising and q = spanning tree cases). 

We should add here a small caveat. We have seen that both for percolation [28] and 
for the general g-state Potts model, when they are not exact, the polynomials Pe(g, t>) 
provide accurate approximations for unsolved problems that improve with an increasing 
number of edges in B. This is the motivation for our conjecture that the polynomials 
converge to the exact critical manifold in the limit of certain infinite bases. However, 
as shown in [28], it is not enough that the number of edges in B tends to infinity. Let 
us say that a generic basis is of size N x M if it contains N unit cells in the vertical 
direction and M unit cells in the horizontal direction. For example, for the kagome 
lattice the unit cell is the six-edge bow tie pattern shown in Fig. |3t, so that the bases of 
Figs. [Tfc and[7Ji are of size 2x3 and 3x2 respectively. In [28], an example was found of 
a lattice for which the lxM basis prediction for increasing M converged to a number 
that, although similar to the numerically known percolation threshold, was ruled out by 
simulations. Thus we should add the condition that in order to produce a sequence of 
algebraic curves of increasing order that converge to the exact critical manifold, both N 
and M must go to infinity simultaneously with a finite aspect ratio 77 = N/M e (0, 00). 
Although it is difficult to probe this at present, it seems at least plausible, if not likely, 
that these infinite basis predictions should also be independent of r\. 

An interesting feature of our findings is that in some sense the computation of the 
polynomials Pb{q, v) acts as a detector of exact solvability for the underlying model. 
To be more precise, it seems that if the polynomial factorises for any choice of basis, 
shedding always the same "small" factor, the zero set of that factor provides (a part of) 
the exact critical manifold. We have seen this mechanism at play in the case of the square 
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lattice, where the small factor is given by ([3]). More generally, for the square lattice with 
checkerboard interactions (see Fig. [2]) the small factor is given by f )22|) . leading us to 
conjecture that this is the exact critical manifold for the checkerboard model. We stress 
that exact results are only expected if the factorisation is systematic, i.e., occurs for 
any B. As an example of the contrary, the 12-edge basis of the kagome lattice led to 
the fortuitous factorisation fl38|) . and we have shown convincingly that neither of the 
factors provide an exact result. We should also notice that the implication does not 
seem to work the other way around: exact solvability does not imply factorisation. This 
is witnessed in particular by the case (q, v) = (3, —1) of the kagome-lattice Potts model, 
which is exactly solvable jlQ], but not a zero of Pb{<1, v) for the finite bases that we have 
studied here. 

Our work hints at several directions for future research. It would obviously be 
worthwhile having a more efficient means of dealing with yet larger bases. This would 
require finding an alternative definition of Ps(q,v) that does not refer to contraction- 
deletion. The study of the critical manifolds for other lattices is another possibility. 
Finally, it is quite possible that Pe(g,f) conceals some graph theoretical applications. 
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Appendix 

Here, we report the graph polynomials Pb(q, v) for the various bases B and graphs G 
considered. 

Square lattice 

The 16-edge basis of Fig. [5b: 

(q - v 2 )(q + Av + v 2 )(q 6 + 12q 5 v + 68gV + 232gV + 4gV 
+ 516gV + 40gV + gV + 736qv 5 + 192gV + 8q 3 v 5 + 576v 6 
+ 504gt> 6 + 52gV + 736v 7 + I92qv 7 + 8gV + 5l6v 8 + AOqv 8 
+ q 2 v 8 + 232v 9 + Aqv 9 + 68v 10 + Uv 11 + v 12 ) (Al) 

The 32-edge basis of Fig. [5fc: 

(q - v 2 )(q + Av + v 2 )(q 14 + 28q 13 v + 384g 1 V + 3424g 1 V 
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+ 4q 12 v 3 + 22240g 1( V + 112q n v 4 + q 12 v 4 + 111744gV 

+ 1536q 10 v 5 + 24g x V + 450016gV + 13632gV + 304g 1( V 

+ 1484032? V + 87520gV + 2688gV + 8q 10 v 7 + 4053840gV 

+ 430560gV + 18520gV + 176gV + q 10 v s + 9198080gV 

+ 1676320gV + 104032gV + 2048gV + 20gV + 17211584gV° 

+ 5249504gV° + 484160gV° + 16288gV° + 256gV° 

+ 25996800gV 1 + 13255488gV 1 + 1867840gV 1 + 98816q 6 v n 

+ 2368gV 4 + 12q 8 v n + 30259968gV 2 + 26606400gV 2 

+ 5923440gV 2 + 479248gV 2 + 17256gV 2 + 208gV 2 + gV 2 

+ 24551424gw 13 + 40761344gV 3 + 15114304gV 3 + 1896032gV 3 

+ 100736gV 3 + 2240gV 3 + 16gV 3 + 10616832w 14 + 43298304gv 

+ 29706752gV 4 + 6060928gV 4 + 484352gV 4 + 16800gV 4 

+ 240gV 4 + 24551424f 15 + 40761344gv 15 + 15114304gV 5 

+ 1896032gV 5 + 100736gV 5 + 2240gV 5 + 16gV 5 + 30259968v : 

+ 26606400gt> 16 + 5923440gV 6 + 479248gV 6 + 17256gV 6 

+ 208gV 6 + gV 6 + 25996800w 17 + 13255488gv 17 + 1867840gV 7 

+ 98816gV 7 + 2368gV 7 + 12gV 7 + 17211584w 18 + 5249504gv 18 

+ 484160gV 8 + 16288gV 8 + 256gV 8 + 9198080w 19 + 1676320gv 

+ 104032gV 9 + 2048gV 9 + 20gV 9 + 4053840*; 20 + 430560gt> 20 

+ 18520gV° + 176gV° + gV° + 1484032w 21 + 87520gv 21 

+ 2688gV 4 + 8gV x + 450016w 22 + 13632gv 22 + 304gV 2 

+ 111744f 23 + 1536gv 23 + 24gV 3 + 22240w 24 + 112gw 24 + gV 4 

+ 3424w 25 + 4qv 25 + 384v 26 + 28w 27 + v 28 ) 



Kagome lattice 

The 24-edge basis of Fig. [7b: 

g 12 + 24q n v + 276g 10 t> 2 + 2016gV + 8g 4 V + 10452gV 
+ 168gV + 40680gV + 1680gV + 122384gV + 10564gV 
+ 28g V + 287760? V + 46440? V + 504? V + 525096gV 
+ 149724gV + 4272gV + 718704gV + 358456gV + 22376gV 
+ 56gV + 673920gV° + 618792gV° + 78864gV° + 864gV° 
+ 331776gw n + 682128gV 1 + 183864gV 1 + 5976gV 4 
+ 270432gi; 12 + 225684gV 2 + 20700gV 2 + 96gV 2 - 248832v 13 

- 95256gv 13 + 10320gV 3 + 264gV 3 - 613008w 14 - 224916gv 14 

- 11592gV 4 - 108gV 4 - 723600w 15 - 143536gv 15 - 3416gV 5 

- 8gV 5 - 550377v 16 - 53844gw 16 - 444g 2 t> 16 - 303288w 17 

- 13512gv 17 - 24gV 7 - 127636w 18 - 2284gt> 18 - 41784v 19 

- 240gv 19 - 10590W 20 - 12gv 20 - 2024w 21 - 276w 22 - 24w 23 - v 2A 

The 36-edge basis of Fig. [7b: 

g 18 + 36g 17 w + 630g 4 V + 7128g lf V + 12g 4 V + 58503g 14 w 4 
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+ 396g 15 w 4 + 370440g 1 V + 6336g 14 t» 5 + 1878670g 12 t> 6 + 65322g 13 v 6 

+ 66q l4 v 6 + 7817940g n w 7 + 486456g 12 w 7 + 1980q 13 v 7 

+ 27122841g 1 V + 2780142g 1 V + 28710g 1 V + 79228584gV 

+ 12642612g 1 V + 267216g 1 V + 220g 1 2u 9 + 195849552g 8 t> 10 

+ 46787952gV° + 1787790g 10 w 10 + 6000g n w 10 + 409851792gV 1 

+ 142829784g 8 w n + 9120348g 9 i; 11 + 78876g 1( V 1 + 722497068gV 2 

+ 361945104gV 2 + 36695373gV 2 + 662258gV 2 + 585g 1( V 2 

+ 1059589512gV 3 + 760782804gV 3 + 118640916gV 3 

+ 3964788gV 3 + 14424gV 3 + 1261019664gV 4 + 1313088660gV 4 

+ 310356054gV 4 + 17868048gV 4 + 169878gV 4 + 54gV 4 

+ 1159847424gV 5 + 1813868856gV 5 + 652644588gV 5 

+ 62170040gV 5 + 1258416gV 5 + 2472gV 5 + 743510016g 2 t> 16 

+ 1892524176gV 6 + 1072172331gV 6 + 167153268gV 6 

+ 6467142gV 6 + 39798gV 6 + 12gV 6 + 251942400gv 17 

+ 1290209472gV 7 + 1270157004gV 7 + 334649052gV 7 

+ 23635824gV 7 + 347148gV 7 + 528gV 7 + 315207936gw 18 

+ 827146404gV 8 + 431966778gV 8 + 57655440gV 8 

+ 1800882gV 8 + 8596gV 8 - 188956800*; 19 - 193284144gt> 19 

+ 116524980gV 9 + 59677656gV 9 + 4411884gV 9 + 51288gV 9 

+ 12gV 9 - 678249936w 20 - 773787816gi; 20 - 180556083gV° 

- 11670198gV° - 299721gV° - 3066gV° - 1195116768w 21 

- 914327676gw 21 - 159944736gV 1 - 8338004gV 1 - 126372g 4 t» 21 

- 336<A; 21 - 1387555272*; 22 - 681761826gw 22 - 75759996gV 2 

- 2362578gV 2 - 17106gV 2 - 12gV 2 - 1197832536i; 23 

- 374476536gw 23 - 25228524gV 3 - 441540gV 3 - 1308gV 3 

- 819069489w 24 - 161325810gw 24 - 6339726gV 4 - 59194gV 4 

- 48gV 4 - 459630612w 25 - 56009760gw 25 - 1218540gV 5 

- 5664gV 5 - 215826246w 26 - 15791340gv 26 - 175608gV 6 

- 360gV 6 - 85595748v 27 - 3593976gw 27 - 17976gV 7 

- 12gV 7 - 28735887w 28 - 648276gw 28 - 1167gV 8 - 8138088w 29 

- 89604gt> 29 - 36gV 9 - 1926908t> 30 - 8946gw 30 - 375648w 31 

- 576gw 31 - 58863w 32 - 18qv 32 - 7140f 33 - 630w 34 - 36w 35 - v 36 (A4) 

The 36-edge basis of Fig. 01: 

g 18 + 36q 17 v + 630q 16 v 2 + 7128g 1 V + 12q 16 v 3 + 58506g 14 w 4 
+ 396g% 4 + 370548g i: V + 6336g 1 V + 1880512g 12 w 6 + 65334g 1 V 
+ 66q 14 v 6 + 7837704^ V + 486876g x V + 1980q 13 v 7 
+ 27272076g 1 V + 2787144g x V + 28728g x V + 80068176gV 
+ 12716088g 1( V + 267888g 1 V + 220g x V + 199480572gV° 
+ 47328210gV° + 1799508g 10 w 10 + 6018q u v 10 + 422114976gV 1 
+ 145763088gV 1 + 9245904gV 1 + 79656g 1 V 1 + 754976880gV 2 
+ 374020200gV 2 + 37614312gV 2 + 676604gV 2 + 612g 1( V 2 
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+ 1126577808g 5 v 13 + 798846408gV 3 + 123476376g 7 v 13 

+ 4118160gV 3 + 15432gV 3 + 1366323552g 4 v 14 + 1404553176gV 4 

+ 329074176g 6 v 14 + 18944598g V 4 + 186240g 8 v 14 + 90g 9 v 14 

+ 1280396160gV 5 + 1977683904g 4 v 15 + 706012560gV 5 

+ 67401848g 6 v 15 + 1413336g 7 v 15 + 3528gV 5 + 834582528g 2 v 16 

+ 2099656080gV 6 + 1181501136g 4 v 16 + 184956942gV 6 

+ 7412292gV 6 + 53976g 7 v 16 + 36gV 6 + 286654464gv 17 

+ 1451444832g 2 v 17 + 1419423912gV 7 + 375629268g 4 v 17 

+ 27409380gV 7 + 457176g 6 v 17 + 1152g 7 v 17 + 358691328gv 18 

+ 931971420g 2 v 18 + 487209582gV 8 + 66732144g 4 v 18 

+ 2296374g 5 v 18 + 15436gV 8 + 6q 7 v 18 - 214990848v 19 

- 21761 1360gv 19 + 128909448g 2 v 19 + 67192668gV 9 

+ 5321868g 4 v 19 + 80928gV 9 + 144g 6 v 19 - 765345024v 20 

- 871306632gv 20 - 208435194g 2 v 20 - 14587674gV° 

- 406632g 4 v 20 - 4050g 5 v 20 - 1334035008v 21 - 1023060564gv 21 

- 184223328g 2 v 21 - 10352660g 3 v 21 - 178764g 4 v 21 - 540g 5 v 21 

- 1528264908v 22 - 755329716gv 22 - 87028104g 2 v 22 - 2969586gV 2 

- 25248g 4 v 22 - 18g 5 v 22 - 1299190104v 23 - 40941 75 12gv 23 

- 28747908g 2 v 23 - 554952g 3 v 23 - 1956gV 3 - 873980685v 24 

- 173604762gv 24 - 7119900g 2 v 24 - 72934g 3 v 24 - 72gV 4 

- 482643540v 25 - 59258256gv 25 - 1341372g 2 v 25 - 6684g 3 v 25 

- 223371438v 26 - 16435896gv 26 - 188820gV 6 - 396g 3 v 26 

- 87528036v 27 - 3687672gv 27 - 18852g 2 v 27 - 12gV 7 

- 29116983v 28 - 657756gv 28 - 1194g 2 v 28 - 8194248v 29 

- 90204gv 29 - 36g 2 v 29 - 1932752v 30 - 8964gv 30 - 376032v 31 

- 576gv 31 - 58875v 32 - 18gv 32 - 7140v 33 - 630v 34 - 36v 35 - v 36 (A5) 



